Time-Domain Measurement of Broadband Coherent Cherenkov Radiation 
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We report on further analysis of coherent microwave Cherenkov impulses emitted via the Askaryan 
mechanism from high-energy electromagnetic showers produced at the Stanford Linear Accelerator 
Center (SLAC). In this report, the time-domain based analysis of the measurements made with a 
broadband (nominally 1-18 GHz) log periodic dipole array antenna is described. The theory of a 
transmit-receive antenna system based on time-dependent effective height operator is summarized 
and applied to fully characterize the measurement antenna system and to reconstruct the electric 
field induced via the Askaryan process. The observed radiation intensity and phase as functions of 
frequency were found to agree with expectations from 0.75-11.5 GHz within experimental errors on 
the normalized electric field magnitude and the relative phase; o~ R \ E \ = 0.039 /iV/MHz/TeV and 
(70 = 17°. This is the first time this agreement has been observed over such a broad bandwidth, 
and the first measurement of the relative phase variation of an Askaryan pulse. The importance of 
validation of the Askaryan mechanism is significant since it is viewed as the most promising way to 
detect cosmogenic neutrino fluxes at E v > 10 15 eV. 

PACS numbers: 84.40.Ba, 95.85.Bh, 95.85.Ry, 41.60.Bq 
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I. INTRODUCTION 

G. Askaryan proposed in 1962 that the charge excess 
in a compact particle shower in a dielectric medium will 
produce a coherent radio Cherenkov emission . Subse- 
quent theoretical work supported this prediction 
The experimental verification came in 2001 [U, with fol- 
low up measurements confirming the frequency and the 
polarization properties of the emitted radiation j(| . 

The interest in the characterization of the Askaryan 
process comes from the idea that it can be used to detect 
ultra-high-energy astrophysical neutrinos {E v > 10 15 eV) 
interacting in radio transparent dielectric media, e.g. , 
ice 0, HM E3, salt @, or Moon's regolith [n], Eg. 
Considering the increasing experimental effort to observe 
such neutrinos, it is of importance to understand the 
properties of coherent Cherenkov radiation by verifica- 
tion of theoretical models. Furthermore, it is of prac- 
tical value to provide the time-domain characterization 
of coherent Cherenkov radiation since the experimental 
triggering and, to large extent, the data analysis will be 
based on time-domain properties of the signal. 

The emission of a coherent radio signal comes from 
the charge asymmetry in particle shower development. 
The asymmetry is due to combined effects of positron an- 
nihilation and Compton scattering on atomic electrons. 
There is ~20% excess of electrons over positrons in a 
particle shower, which moves as a compact bunch, a few 
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cm wide and ~1 cm thick, at a velocity above the speed 
of light in the medium. The frequency dependence of 
Cherenkov radiation emitted is dP cx vdv. For radiation 
with wavelength A ^> I, where I is the length scale of 
the particle bunch, the radiated electric field will add co- 
herently and thus be proportional to the shower energy. 
A radio signal emitted by a particle shower in a dielec- 
tric material such as ice or salt is coherent up to a few 
GHz, is linearly polarized, and lasts less than a nanosec- 
ond. A shower with energy of 10 19 eV interacting in 
the ice will produce a radio pulse with peak strength of 
~ 1 mV/m/MHz at the distance of 1 km. 

This report describes the analysis of measurements of 
of Askaryan pulses recorded with a log periodic dipole 
array (LPDA) antenna in an experiment (SLAC T460) 
performed at the Final Focus Test Beam (FFTB) facility 
at SLAC in June 2002 |6J. The data taking is described 
first, followed by the discussion of a time-domain based 
description of the antenna system response to transient 
radiation, and the description of system calibration. In 
the last two sections, the reconstruction of the electric 
field induced by the particle shower is presented and com- 
pared to the theoretical work. 



II. MEASUREMENTS 

The experimental setup is illustrated in Fig. ^ which 
shows top and side views of the salt-block target, re- 
ceiving antenna locations, and the incoming beam. The 
particle showers in the salt target were induced by 
bremsstrahlung gamma-ray photon bunches emanating 
from an aluminum radiator plate placed in the path of 
the electron beam. The electron beam itself was deflected 
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FIG. 1: Geometry of salt-block target and receiving antennas in SLAC T460 experiment, from Ref. 



to pass below the salt target. The electrons in the beam 
were accelerated to ^28.5 GeV and the energy of sec- 
ondary bremsstrahlung photons followed E~ x distribu- 
tion up to ~28.5 GeV. The total energy of photons in 
each bunch varied with the electron beam intensity and 
the thickness of the radiator plate. 

The salt target was constructed from salt bricks, with 
a total mass of about 4 metric tons. The detailed descrip- 
tion of the target construction can be found in Gorham 
et al. [||. The radio pulses arising from photon induced 
particle showers were collected by bowtie antenna array 
embedded in the salt and by a C/X-band horn antenna 
and an LPDA antenna outside the salt. The bowtie an- 
tennas were located in a plane parallel and about 35 cm 
above the shower axis. The horn and LPDA antennas 
were placed in the plane of the shower, with LPDA point- 
ing approximately at the expected shower maximum, as 
indicated in the Fig. ^ The side of salt target facing 
the external antennas was angled at 10° in order to allow 
transmission of Cherenkov radiation which would oth- 
erwise have been totally internally reflected. The an- 
tennas were read out using digital oscilloscopes, and to 



improve signal-to- noise ratio (SNR), the recorded wave- 
forms were averaged over many pulses, using an ultra- 
stable microwave transition-radiation trigger from an up- 
stream location. 

The majority of data collected in the experiment was 
analyzed and published in Ref. ||. In order to fur- 
ther test theoretical models of the coherent Cherenkov 
radiation, this work will concentrate on the analysis of 
the data collected with the LPDA in two runs, 35 and 
109, in which no microwave filters were placed between 
the LPDA and the oscilloscope, 1 providing the broad- 
est bandwidth data. The total photon bunch energies in 
these runs were estimated to be 0.48 and 1.1 EeV, re- 
spectively, with an uncertainty on the order of 20% due 
to the beam intensity and the electron bunch energy dis- 
tribution fluctuations. 

The LPDA is an antenna constructed by arranging 
multiple dipolcs in a geometry that provides high signal 



1 Although, a 20 dB attenuator was used in run 109 to restrict the 
peak voltage to the dynamic range of the oscilloscope. 



3 




-1 - 



b 2 4 6 8 10 12 14 16 18 20 
time, ns 




2 4 6 8 10 12 14 16 18 20 



frequency, GHz 




frequency, GHz 

FIG. 2: (a) Raw voltage recorded with the LPDA in run 109, 
and its (b) magnitude and (c) group delay as functions of 
frequency. 

gain over a large frequency bandwidth 01 • The LPDA 
antenna used was Electro-Metrics model EM-6952 with 
nominal bandwidth from 1-18 GHz. The antenna was 
connected by two pieces of 75-foot heliax cable, Andrew 
LDF4-50A, and by three pieces of 12-inch semi-rigid 
Haverhill cable, to a CSA8000 sampling oscilloscope with 
20 GHz bandwidth and up to 1000 GSa/s sampling rate. 

The raw signal recorded in run 109 is shown in Fig. |2 
along with its amplitude and group delay, defined as 
Tg{y) — (27r)~ 1 , where 4>{v) is the relative phase 
of the voltage frequency component, <fi = tot. It can be 
seen that the signal is present from ~0.75 GHz (where the 
LPDA loses sensitivity) to ^7.5 GHz, where the signal 
strength drops down to the intrinsic oscilloscope noise 
level. Run 35 signal looks identical except for the in- 
creased bandwidth (up to ~11.5 GHz) since no attenua- 
tor was used. Fig. shows spectrograms of the recorded 
signals, with the curvature indicating variations in the 
delay of voltage frequency components due to the LPDA 
and the transmission line. 
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FIG. 3: (Color online) Spectrograms of signals recorded in 
runs 109 and 35 (top and bottom). The color indicates the 
logarithm of the square root of the signal power per unit fre- 
quency. 



III. SYSTEM CHARACTERIZATION 

In order to reconstruct the electromagnetic radiation 
pulse incident onto the LPDA (or any other antenna), 
the signal distorting effects of the antenna and the trans- 
mission line have to be corrected. The procedure out- 
lined here is general and applicable to reconstruction of 
any observation of impulsive radio signals. The antenna 
equations and the required RF parameters are described 
in the appendix B 



A. Antenna system calibration 

The measurement of the LPDA effective height was 
performed in an anechoic chamber at University of 
Hawaii by a reciprocal S21 method. Two identical anten- 
nas were mounted about 60-in apart facing each other, 
which ensured they were in each other's far-field re- 
gion. For an LPDA, the far-field requirement reduces 
to d > A/4, where d is the distance from the point at 
the LPDA where radiation of wavelength A preferentially 
couples, i.e., the phase center. The transmitting antenna 
was stimulated by a 200-mV step impulse generated by 
an HP 54121A logging head. The received signal was 
amplified by an Agilent 83017A broadband amplifier and 
recorded by an HP 54120B digitizing oscilloscope with 
a 20 GHz bandwidth at 100 GSa/s. The presence of an 
amplifier introduces a transfer function to the system, so 
that the expression for the recorded signal V rec , combin- 
ing Eas. IA.10l axid lA.lll is 

Vrec(t) = T^—H amp {t) o h N {t) o h N {t) o dV ^f \ (1) 
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where r is the antenna separation, c is the speed of light, 
V src is the voltage stimulating the transmitting antenna, 
hjsf is the normalized effective height vector of the an- 
tennas, and H amp is the transfer function correction ac- 
counting for cables and the amplifier. The scalar product 
(see Eq. lA.lOj) is omitted since the effective height vec- 
tors were parallel. The quantity H amp was measured by 
the same setup, but excluding the antennas from the cir- 
cuit. Similarly, the heliax cable used in the SLAC mea- 
surement was stimulated on one side by a 200-mV step 
and the resulting pulse was recorded at the other end 
with the digitizing oscilloscope. The semi-rigid Haver- 
hill cable was unavailable for time-domain transfer func- 
tion calibration, so only the attenuation was measured 
as a function of frequency with a network analyzer. Its 
phase response was ignored, but due to the relatively 
short length of that cable this omission will have a very 
small effect on the final result. 



B. Mathematical operations 

In order to calculate the LPDA effective height from 
calibration measurements, it is conceptually easiest to 
take a Fourier transform of Eq. ^ an d re-arrange it so 
that 



h N (u) 



V rec {v) 



iv H amp (v)V src {v) 



(2) 



where f{v) = J f(t)e 27Tll,t dt. Also, from Eq. |A~TT1 
Hamp(v) = Vout(y)jVin{v). The same expression is also 
applied to calculate the transfer function of heliax cable. 
All operations, including the square root, are performed 
using complex quantities. Special care should be taken 
when taking a complex square root in order to properly 
account for phase wrapping. 

While straightforward, this frequency-domain based 
approach will give poor time-domain results if the fre- 
quency bandwidth of a measured signal is substantially 
narrower than the Nyquist limit of the sampling oscil- 
loscope. The time-domain result of a complex division 
operation in the frequency-domain will be dominated by 
frequency components where the signal is absent, i.e., 
where a division of two small noise values can produce 
an arbitrarily large artificial result. Digital signal pro- 
cessing tools can be applied to filter out such artificial 
"out-of-band" noise, but they will in all cases produce 
an unacceptable level of signal distortion. In order to 
circumvent this issue, a time-domain deconvolution algo- 
rithm with noise reduction can be applied in place of ev- 
ery complex division o per ation. In this work, the Wiener 
algorithm was chosen [Tj|, and the level of the noise-to- 
signal power ratio required by the algorithm was chosen 
such to maximize out-of-band noise rejection while pre- 
serving the fidelity of the in-band signal (see Fig. |Sp) . 

The resulting time-domain response of the LPDA-cable 
system as used in the experiment can be expressed as 
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FIG. 4: (a) Time-domain response, and (b) attenuation and 
(c) group delay as functions of frequency for the LPDA-cable 
system used in the measurements. The phase coherence of 
calibration pulses extends to ~13 GHz, which exceeds the 
usable bandwidth of experimental data. The feature visible 
in group delay plot at ~8.9 GHz is due to reflections at the 
amplifier input and it was not possible to calibrate it out. 
However, it does not create any appreciable effect in the final 
result. 



h sys (t) = H cab ie(t) o h N , LPDA (t) and is shown in Fig.gl 



IV. ANALYSIS 

The voltages recorded in the experiment can be ex- 
pressed as (see Eqs. IA.6HA.8(l 



Vrec(t) 



h sys (t) * E ant (t), 



(3) 



where E ant is the electric field at the antenna due to the 
particle shower. During the measurements, the LPDA 
was aligned such that the polarization of electric field 
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FIG. 5: The electric fields at the antenna in two runs (a), 
and the magnitude (b) and the group delay (c) as functions 
of frequency (for run 109 only). In (b), the dashed line shows 
the magnitude obtained by straight frequency-domain divi- 
sion. The signal is dominated by artificial noise at frequen- 
cies outside of the signal bandwidth (0.75-7.5 GHz) . It is clear 
that the use of the Wiener deconvolution algorithm is justi- 
fied since it preserves the signal "in-band" while suppressing 
the "out-of-band" noise. 



was parallel to the effective height vector, so that the 
spatial scalar product drops out of the equation. Thus, 
in the frequency domain, the electric field at the antenna 
is given by 



E an t{v) 



Zv V rec (v) 



Z L h s 



(4) 



with Zl — 50fi. The resulting electric field is shown in 
Fig. El 

Before the measured electric fields can be compared 
with the theoretical expectation, several corrections have 
to be applied. The theory was derived for a pulse de- 
tected at a very large distance from a shower initiated by 
a single particle, i.e., where the shower can be considered 



a point source emitter and the radiation emitted is coher- 
ent over the full length of the shower. Additionally, the 
detection is also assumed to be in the same medium as the 
emission. In the experimental setup this is not the case. 
Aside from accounting for the electric field divergence 
as it crosses from salt to air, the portion of the particle 
shower generating a coherent pulse at the antenna scales 
with the frequency of observation. Finally, the number 
of particles at the shower maximum, N max , due to the 
superposition of many low-GeV electromagnetic showers 
will be larger than N max due to a shower initiated by 
a single particle of the equivalent energy. Since the in- 
tensity of the radiation emitted is proportional to N max , 
this difference has to be taken into account. 



A. Coherence zone correction 

The standard radiation coherence requirement is that 
phases of rays emerging from two ends of a coherence 
zone do not differ by more than one cycle (2tt) at any 
given frequency. The geometry of this requirement in 
the present case is sketched in Fig.|SJi. Thus, the zone of 
coherence about a point on the shower axis that the an- 
tenna is pointed at (zc), can be defined for z's satisfying 
the condition 



R'i{z- v )-R'i(zc; V ) 



zc 



2tt 

T 



< 7T, 



(5) 

where zq — 0.53 m, n = 2.44 is the index of refraction 
of salt, radiation wavelength A — — , and R"(z; v) is the 
apparent distance from a point on the shower axis to the 
phase center of the antenna for the given frequency as 
if the entire path was in salt, i.e., the number of actual 
phase cycles, N, is preserved and R" = NX. Eq. [5] is 
simply a restatement of the difference in phase factors of 
an electric field due to a moving charge without making 
a Fraunhofer zone approximation |15| . For the phase 
center of the LPDA, a simple model is assumed where 
the phase center moves linearly along the antenna axis, 
such that the 1-GHz phase center is at the feed-point 
and the 20-GHz phase center is at the tip of the antenna, 
about 26 cm from the feed-point. While not exact, this 
model is sufficiently accurate not to contribute strongly 
to systematic errors. As an example, the coherence zones 
at 1 GHz and 5 GHz were found to be z = [0, 1.05] m 
and z = [0.45, 0.77] m, respectively. 

With the coherence zone defined, a correction to the 
electric field at the antenna with respect to the electric 
field at infinity can be made. From Refs. 0, 0> 0; the 
electric field due to coherent Cherenkov radiation at a 
given angle is proportional to the integral over the current 
distribution in the particle shower, J(z) — c q(z), and 
to geometric factors relative to emission and observation 
points, 



E{v,e) 



ik(R cos 9-\-z/n) 



q(z)dz, 



(0) 
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FIG. 6: a) A sketch of the experimental geometry, including 
the coherence zone (shaded) and the charge excess distribu- 
tion as measured by the embedded bowtie antennas |£|. The 
dashed line, labeled virtual shower image, indicates the non- 
trivial path along which the shower appears to develop due 
to the salt /air interface, as seen by an observer located at the 
LPDA feedpoint. b) The correction factor magnitude and c) 
phase as functions of frequency (see text). The dashed line in 
b) is the ad-hoc correction factor used in Ref. @. 



Cm t g 



sin 8c 



dz 



Jo q{z)dz 



(7) 



where [z\{v), defines the coherence zone at a given 

frequency, T = T(z;v) is the transmission coefficient 
through the salt /polyethylene/air interface, G = G{z\v) 
is the relative gain of the LPDA away from the antenna 
boresight, and = 6(z\v). The last two terms of the 
integrand in the numerator account for the drop in the 
electric field magnitude away from Cherenkov angle, with 
AO = 4.65° [i^s] 0. 

The transmission coefficient is given by the Fresnel 
equation for the E-field parallel to the plane of inci- 
dence [ljl, 



rp M. 

~ En 



2ncosi 



n cos r + n' cos i 



where i and r are incident and refracted angles related 
by Snell's law, and n and n' are indices of refraction in 
the two media. In calculating the transmission coeffi- 
cient, the presence of the 2.5 cm thick polyethylene sheet 
{n po iy = 2.25) is ignored at frequencies below 2.5 GHz 
(half-wavelength thickness) and its effect is gradually 
turned on to full strength above 10 GHz (two-wavelength 
thickness) . The relative gain of the LPDA as a function 
of frequency was not measured, but was modeled based 
on the manufacturer's specification of 70° average beam- 
width in the antenna E-plane 3 and the standard LPDA 
expectation of first null at twice the beam-width, to be 

G(V>) ~ y cos(|g^), with %p being the angle relative to 

the antenna boresight. 

When presenting the preliminary results of this anal- 
ysis in Ref. Q, it has been argued that the coherence 
zone extends over the full shower at frequencies below 
800 MHz, and that at higher frequencies the near-field 
correction is proportional to v. In Fig. 03, it can be seen 
that this simplification was reasonable as it produced an 
adequate near-field correction. 



where k = ^ , R is the distance from an emission point to 
the observation point, 9 is the angle of emission, and the 
integral is taken over the full length of the shower. Taking 
into consideration the geometry of the experiment, 2 a 
near-field correction factor relative to observation made 



B. Field divergence at the salt /air interface 

The 1/R attenuation in the electric field has to be mod- 
ified for the field divergence at the salt/air interface. The 
field divergence factor for a spherical wave incident on a 
planar interface can be calculated by considering an area 
element dA of the surface, 

R 2 s dn R' 2 dQ' 
dA = —2 — = —5 , (8 

cos i cos r 



2 For ease of calculation and presentation, all geometry dependent 
quantities will be parametrized in terms of z, the distance along 
the shower axis, and v, the frequency. 



3 The plane parallel to h. 
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where d£l is the solid angle element subtending dA as 
seen from the emission point, d£l' is the solid angle el- 
ement subtending field lines exiting through dA, R s is 
the distance from the emission point to the salt/air in- 
terface, R' s is the distance from the virtual emission point 
as seen from the outside, and i and r are incident and 
refracted angles as determined by Snell's law. The field 
divergence factor is given by the ratio of true and virtual 
distances 0, 



R s n cos i 
R{ cosr 



(9) 



The effective distance from an emission point to the 
LPDA is then given by, 

R'{z;v)=R a {z;v) X±J. +Rb ^ z . v ) i 

cosr(z; v) 

where R a is the distance from the salt / air interface to the 
LPDA. The reader can easily convince oneself that the 
presence of the polyethylene sheet cancels out when con- 
sidering the field divergence in the experimental setup. 
The mean LPDA to shower distance at any given fre- 
quency is then 



(R'{v)) = 



N„ 



SllR\z-v)dz 
f dz • 

J z-i 



correction 



(ii) 



The energy of photons initiating particle showers in the 
SLAC experiment was from -1-28500 MeV, with al/£ 
bremsstrahlung distribution. The theoretical work on co- 
herent Cherenkov radiation is based on simulations of 
electromagnetic showers initiated by primaries with en- 
ergies of 1 TeV and above 0,0 . While N max , the number 
of particles at the shower maximum, scales nearly linearly 
with the energy of the primary at E p > 10 TeV, this is 
not the case for the lower energy primaries. The EGS4 
simulation shows that a bremsstrahlung photon bunch 
with the total energy of 20.7 TeV hitting the salt tar- 
get will produce N max w 55000. However, the 20.7-TeV 
electron primary, interacting in the salt, gives rise to a 
shower whose peak is deeper and broader, and hence has 
only —19500 particles at the shower maximum [lrji].Thus. 
CNmax = 2.82 can be defined as a correction factor to be 
used when comparing to theoretical results. 

Putting all these corrections together, the measured 
electric field can be compared to the one which would 
have been expected from an equivalent shower observed 
in salt in the far-field region at Cherenkov angle, 



R\E(v,0 c )\ 



{R'{v))E ant (v) 

CnfCN III (!.r 



(12) 



The corrections will be applied only in the frequency 
range where the coherent signals are present in order not 
to amplify the thermal noise background. 
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FIG. 7: (a) Electric field due to Askaryan pulses from two 
runs renormalized to 1 TeV primary energy; (b) electric field 
magnitudes as functions of frequency compared to the ex- 
pectations; (c) electric field phases as functions of frequency, 
calculated relative to the zero-crossing point, compared to the 
expectation for an Askaryan pulse in salt (solid line) [3. The 
dashed line places the theoretical phase expectation in the 
fourth quadrant, see text. 



V. DISCUSSION 

The far-field corrected pulses, and their magnitudes 
and phases, are shown in Fig. [7| The magnitude and 
phase information have been binned in 250 MHz bins 
for clarity, and error bars indicate the RMS variation 
within each bin. The expectation for the magnitude of 
the electric field is given by 0, @ 



R\E\ = A f d 



E p 




V 






ITeV 








1 + 



1 



U-44 



(13) 



where A = 2.53 x 10~ 7 V/MHz, f d = 0.52, v = 
1.15 GHz, and v x = 2.86 GHz. The standard Fourier 
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transform 4 of the measured electric field was multiplied 
by a factor of 2 in order to agree with the definition of 
Eq. 8 in Ref. adopted in Refs. 0,0. 

The only published electric field phase expectation has 
been calculated for an Askaryan pulse in ice 2] , but there 
is no reason to expect any significant deviation in the 
case of salt. However, there are two ambiguities that 
have to be resolved before a meaningful comparison of 
measured and expected phases can be made. The first 
is the actual quadrant of the phase. In the early theo- 
retical calculations, the signs of the real and imaginary 
parts of the Fourier transform of the electric field were 
not kept |17| . and thus the complete phase information 
was unavailable. Newer calculations indicate that the 
phase is expected to lie in the fourth quadrant |17|. and 
the dashed line in Fig. 0: indicates such a choice, which 
will be considered as correct in the remainder of this dis- 
cussion. The second ambiguity is in the choice of the 
reference time about which to calculate the phase angle. 
No convention has been put forward, but one could ex- 
pect that a natural choice would be the time of arrival 
of radiation from the shower maximum in the direction 
of the Cherenkov angle. Unfortunately, in the present 
experiment the relative oscilloscope trigger time with re- 
spect to the beam pulse was not recorded. The best guess 
would be to take the zero-crossing time in the middle 
of the pulse as was done in Fig. 0:. A choice of a dif- 
ferent reference time introduces an artificial phase slew, 
which can obscure the real physical features related to 
the phase, like the current distribution in the shower. 

A few features of the recorded pulses should be noted. 
The signals were recorded from 0.75-11.5 GHz in run 
35 and 0.75-7.5 GHz in run 109. This is indicated by 
the loss of both the phase coherence and the SNR above 
these frequencies due to attenuation in the antenna sys- 
tem, which increases with the frequency As noted before, 
a 20 dB attenuator was present in run 109, causing the 
additional bandwidth loss. A clear way to present an im- 
pulsive, coherent, broadband signal expected from coher- 
ent Cherenkov radiation is with a spectrogram (Fig. 0. 
The pulse arrival is seen at zero time, with some low level 
reflections trickling in at the later times. The effect of 
deconvolving system response, Eq. 01 and applying cor- 
rections, Eq. El can be seen by comparing to Fig. [3] 

Although deviations from theoretical expectations in 
Figs. Et> an d Et appear to be systematic since they 
are correlated between the runs, they are actually due 
to noise in the calibration data. This can be clearly 
seen in Fig. [§] which histograms bin-by-bin differences 
between expectation and measurements over the signal 
bandwidth. The full frequency resolution of 50 MHz 
per bin was used. The deviations from expectation are 
Gaussian in nature, implying random noise. The off- 
sets in the means of difference in magnitude can be at- 



log 10 (V/m/MHz) 




FIG. 8: (Color online) Spectrograms of electric fields due to 
Askaryan pulse recorded in runs 109 and 35 (top and bottom) , 
where the color indicates the logarithm of the square root of 
electric field power per unit frequency. The reconstructed 
impulses are coherent and unresolved at all frequencies. 



Run 35 



Run 109 




E(tu) = jr oo E(t)e^dt. 



FIG. 9: Distributions of bin-by-bin differences in expected 
and measured magnitudes and phases of electric fields, over 
signal band widths noted in the text. The magnitudes have 
been renormalized to 1 TeV shower energy. The fits are Gaus- 
sian. 



tributed to systematic errors in the energy calibration 
of the runs, while the offset in the means of difference 
in phase, can be attributed to the choice of phase refer- 
ence time. Changing the reference time by 8(1) ps for 
run 35(109) shifts the mean to 0° while preserving the 
width of the distribution. Considering that the digitiza- 
tion time resolution is 10(20) ps, this indicates a shift of 
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less than one timing bin in both cases. Also, one could 
chose to use the predicted intensity of Cherenkov radi- 
ation in order to calibrate the total energy of showers 
in two runs. Minimizing the mean of the difference in 
magnitude, the shower energies are found to be 0.38 EeV 
and 1.36 EeV, respectively. Combining the distributions 
from two runs, the fitted deviations are A(R\E\/E P ) = 
-0.005 ± 0.039(stat) ± 0.045(sys) ^V/MHz/TeV and 
A<j) = -9° ± 17°(stat) ± 22°(sys). Thus, it can be con- 
cluded that the theoretical expectations for magnitude 
and phase of the electric field 0, El due to coherent 
Cherenkov radiation in salt in the frequency range from 
0.75-11.5 GHz have been confirmed within experimental 
uncertainty. 

In this work, the procedure for the time-domain based 
analysis of antenna measurements and calibrations has 
been described and successfully applied to measurements 
of coherent Cherenkov radiation made at the FFTB facil- 
ity at SLAC. The time-domain signal analysis preserves 
both magnitude and phase information of the original 
pulse, allowing for more accurate validation of theoreti- 
cal models. The most precise validation of electric field 
intensity and the first validation of electric field phase of 
the theoretical model has been performed. 



APPENDIX: ANTENNA EQUATIONS 

The voltage at the antenna port delivered to the trans- 
mission line of characteristic impedance Zl by the an- 
tenna with impedance Z a R X can be defined as fill Il9). 



with 



V L {t) = T Rx h e ,Bx(Jt)*E{t), 



(A.l) 



with 



2Zt 



TRx 



Z, 



a,Rx 



where * is the operator combining temporal convolution 
and spatial scalar product, 5 E is the incident electric 
field and h e ^R X is the time-dependent effective height of 
the antenna for reception. The time-dependent effective 
height is the time-domain representation of the antenna 
complex impedance, and is equivalent to the antenna re- 
sponse to a delta-like impulse. 

At this point it is useful to also write down the equa- 
tion for the electric field at some distance from the an- 
tenna transmitting an impulsive signal driven through 
the transmission line of characteristic impedance Zl by 
a voltage source V src , 

E(t) = -r^—frx t Tx h e . Tx {t) o V„ e (t), (A.2) 
iirrc 



ITx 



TTx 



Zq 

Z a ,Tx ' 

2Z at Tx 
Zl, + Z n _l 



where o is the convolution operator, Z a j<x is the trans- 
mitting antenna impedance, Zq is the impedance of free 
space (377S1) , and h e ^Tx is the time-dependent effective 
height of the antenna for transmission. The voltage V src 
as used here is the voltage that would have been read by 
an oscilloscope if it were connected to the voltage source 
instead of an antenna and matched to the transmission 
line. 6 Eq. IA.2I can be rewritten using the relation be- 
tween transmission and reception effective heights of an 
antenna derived from self- reciprocity arguments |l8j|. 7 



h e ,Tx(t) = 2d t h eiRx (t). 



(A.3) 



Substituting Eq. IA.3I into Eq. IA.2I and noting that time 
derivation and convolution commute, the transmission 
equation can be written as, 



E(t) = 



' hxTrxK^O^lM. (A.4) 



2irrc 



dt 



Both antenna transmission and reception are now defined 
using the same effective height quantity. 

The form of the preceding equations assumes that both 
the antenna and the transmission line have purely resis- 
tive impedance. The equations can be simplified and the 
complex antenna impedance can be reintroduced by ab- 
sorbing the voltage transmission coefficient into the defi- 
nition of effective height and rcnormalizing voltages and 
electric fields as suggested by Farr and Baum [2(j ■ Using 
the following variable substitutions, 



2\/ Z Zl 




e,Rx 



(t), 



(A.5) 
(A.6) 
(A.7) 



where Z ai V, and E stand for both transmission and 
reception cases, the antenna equations become 



V N<L {t) = h N (t)*E N (t), 



E N (t) = 



1 



2nrc 



h N {t)- 



dV~N,src{t) 

dt 



(A.8) 
(A.9) 



f(t)*g(t)sjf(t)-g(t-t')dt'. 



6 Thus, V sr c is the voltage delivered to the antenna, not the open 
circuit voltage of the source. 

7 The sa me conclusions based on self-reciprocity are derived by 
Baum lltil . Although it would appear that his result, h e ^Tx = 

Rx (using th e no tation adopted in this paper) is in disagree- 
ment with Ref. |l8j. that is not the case since he uses slightly 
different definitions for various antenna system quantities. 
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Finally, if two identical antennas separated by distance r 
are used to transmit and receive, and transmission lines 
with the same impedance are used at both ends, the volt- 
age observed on the port of the receiving antenna is given 
by 

V L {t) = ^h N {t)*h N {t)o^^. (A.10) 

The transmission line complex impedance can be han- 
dled by introducing a transfer function correction, H, 
which accounts for phase and amplitude distortions of a 
broadband signal traveling through the transmission line. 
In the time-domain based description, this correction is 
expressed as 

V out (t) = H(t)oV in (t). (A. 11) 

All of the antenna equations above need to include 



transmission line transfer function corrections at appro- 
priate places, depending on the actual system setup. 
One should keep in mind that convolution operators 
commute, so that contributions from several transmis- 
sion line segments (even ones at the opposite transmis- 
sion/reception ends of the system) can be combined into 
a single transfer function. 
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